Frohlich mass in GaAs-based structures 
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The Frohlich interaction is one of the main electron-phonon intrinsic interactions in polar materials 
originating from the coupling of one itinerant electron with the macroscopic electric field generated 
by any longitudinal optical (LO) phonon. Infra-red magneto-absorption measurements of doped 
GaAs quantum wells structures have been carried out in order to test the concept of Frohlich 
interaction and polaron mass in such systems. These new experimental results lead to question the 
validity of this concept in a real system. 
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Based on the Frohlich interaction, polaronic effects 
were early studied by Lee and Pines 1] and also by Feyn- 
man Q on theoretical grounds: simple consequences were 
derived leading to the concept of the polaron mass, and 
reproduced in many text books This conceptual in- 
teraction also attracted much attention, on the experi- 
mental side, in either three or quasi two-dimensional (3D 
or Q2D) semiconductors Q which were model systems to 
study the effect. For the simplest case of a Q2D semicon- 
ductor, the single-particle energy spectrum of electrons 
of mass m* , in a perpendicular magnetic held B± , con- 
sists of equally spaced (hcu c j_ = eB±/m*) energy states 
known as Landau levels (LL) labelled with an integer in- 
dex n going from to infinity: one direct consequence 
of the polaronic concept is the prediction of the reso- 
nant magneto-polaron coupling (RMPC) which antici- 
pates a non linear field variation of the cyclotron res- 
onance (CR) energy, Hlo c , when hto c ~ Huilo, with a 
pinning of the CR above Tilolo- Such experimental ev- 
idence is quite difficult to obtain because the range of 
energies concerned coincides with that of the strong ab- 
sorption of the phonon band (Restrahlcn) of the material, 
therefore obscuring the data. The theory itself raises also 
fundamental questions because, in real materials, one al- 
ways has many electrons which are known to oscillate 
in phase as plasmons developing in turn a macroscopic 
electric field which couples to that of the LO phonon, 
generating hybrid plasmon-phonon modes. What really 
happens? Does the polaronic concept remains valid? To 
our knowledge there is no quantum mechanical treatment 
of the hybrid modes and therefore it is difficult to answer 
these questions on theoretical grounds. 

Experiments performed in suitable conditions are how- 
ever able to bring new light in that matter. Overcom- 
ing the problem of the Restrahlen band is indeed tech- 
nically very difficult in 3D semiconductors. In that re- 
spect the Q2D structures, grown by molecular beam epi- 
taxy (MBE), are more appealing because, besides the 
fact that, with the remote doping technique, carriers can 
remain quasi-free at the lowest temperatures, one can 
remove, by selective chemical etching, the whole MBE 



structure from the substrate (lift-off process). It is then 
possible to deposit it on an infra-red transparent sub- 
strate, like Si for instance, which is wedged with an angle 
of few degrees to avoid interference effects. This struc- 
ture is composed of a single doped GaAs quantum well 
(QW) of width L and areal density n$ sandwiched be- 
tween two GaAs-AlAs superlattices with the Si-n type 
doping performed symmetrically on both sides of the QW 
[j| . This enables to grow structures of quasi 2D electron 
gas (2DEG) combining high density and high mobility. 
We already reported results on such a sample with 
L = 10 nm,ns = 12 x 10 n cm -2 Q. The same process 
has been applied to new samples with similar structures 
but with L = 13 nm and various lower doping levels. 
Magneto infrared transmission measurements are done 
in an absolute way for any fixed magnetic field B: this is 
achieved with a rotating sample holder allowing to mea- 
sure reference spectra in the same conditions on a par- 
ent wedged piece of Si substrate. With the k vector of 
the incoming light parallel to ~3, the growth axis of the 
sample was maintained either parallel to this direction 
(perpendicular Faraday (PF) configuration) or tilted by 
an angle 9 with respect to it (tilted Faraday (TF) con- 
figuration). The light was guided through an over-sized 
guide pipe ending by a cone to minimize the divergence 
of the beam. The interpretation of the experimental re- 
sults is done by comparing them with a simulation of the 
multi-layer transmission by an appropriate model. This 
is essential because, in the frequency range of interest, the 
spectra can be distorted by dielectric interference effects 
independent of any specific electron-phonon interaction. 

The previous study [(J reported on PF configuration 
measurements where no sign of RMPC was observed. 
This was not necessary a proof that the effect was not 
existing because it could have been screened by the high 
value of ns : indeed theoretical simulations includ- 
ing the RMPC, of data obtained with such high density 
QW was still showing a small interaction. New exper- 
iments with a lower ns, for which this simulation pro- 
vided clear effects, were necessary to check the concept 
of the RMPC. The doping level of the new lift-off samples 
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FIG. 1: Experimental transmission spectra (empty dots) of 
sample 1200. a) absolute transmission for B = 0T and 8 = 
60°; b) and c) absolute transmission and relative transmission 
respectively for B±_ = 20.7 T and 8 = 20°. The continuous 
red curves are the fit of data with the model described in the 
text. 



ranges between ns = 7.0 x 10 11 cto -2 (sample 1200) and 
ns = 9.0 x 10 11 cm~ 2 (sample 1211) with respective mobil- 
ities of /idc — 260m 2 /V/sec and /«dc — 280m 2 /V/sec 
0. All experimental results reported here have been 
obtained at a fixed temperature of 1.8 K. We will con- 
centrate on the results obtained for B± around 20 T for 
which the CR span the Restrahlen band of GaAs. For 
both samples, in that range of fields, the filling factor 
of the structure v — nsh/(eB) remains essentially lower 
than 2 and, in that case, non-parabolicity (NP) effects 
can be ignored. 

Typical spectra, for sample 1200, are displayed in 
Fig.l. The absolute transmission (AT) at B = 0T and 
6 = 60° is shown in Fig. la. The strong absorption 
features are those related to the TO phonons of GaAs 
and AlAs whereas the weak ones are those related to 
the LO modes which are known to become active in ab- 
sorption for thin slabs in the TF configuration. They 
clearly originate from the layers of the superlattices sur- 
rounding the doped QW. They are slab modes and give 
a lower bound of the energy of LO modes in the whole 
structure. In Fig. lb, the absolute transmission of the 
sample in a TF configuration is displayed for 9 = 20° 
and B± = 20. 7T. In practice with the available mag- 
netic fields, we are led to work with angles where the 
LO phonons are less visible but they are still clearly 
present as indicated by the vertical dashed lines. Two 
CR transitions, labelled CR1 and CR2 (vertical dotted 
lines), are now present, the CR2 feature being, for that 
field, at a lower energy than hu>Lo(G&As). In Fig. lc we 
show the relative transmission (RT) obtained by divid- 
ing the AT spectrum of Fig. lb by the AT spectrum at 
B = 0T for the same angle. The RT spectra are a combi- 
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FIG. 2: Observed cyclotron energies for sample 1200 (a) and 
sample 1211 (b) for different tilt angles . The crosses and 
black down triangles are the measured LO and TO energies 
respectively and the empty squares the cyclotron energies in 
the PF configuration. Other points are the observed CR en- 
ergies in the TF configuration: (a) full dots : 8 = 14°, empty 
dots : 8 = 20°, stars : 8 = 25° and (b) full dots : 8 = 15°, 
empty dots : 8 — 26.5°. Results in (b) for B± lower than 
18 T correspond to filling factors > 2 when non parabolicity 
effects are apparent. 



nation of four experimental ones and therefore the noise 
is more important. However the only remaining features 
are those related to the contribution of the 2DEG. The 
energies of the CR1 and CR2 structures for each sam- 
ple, traced as a function of B±, are displayed in Fig. 2 
for different angles. Whereas in the PF configuration 
(empty black squares), no sign of interaction is observed, 
in the TF configuration a pronounced anticrossing of the 
CR1 and CR2 components occurs, which increases with 
(9, one of the components remaining pinned at an energy 
of 35.6 ± 0.05 meV for the sample 1200 and 35.5 ± 0.05 
meV for the sample 1211, both energies being lower than 
?ia;Lo(GaAs). Other samples with different densities of 
electrons (ns = 5.8 x 10 n cm~ 2 and ns = 8.8 x 10 n c77i~ 2 
and the same L of 13 nm or with a smaller L of 10 nm 
and ns — 12 x 10 u cto~ 2 , have also been investigated and 
demonstrate exactly the same qualitative features. 

To interpret these results, one has to calculate the 
transmission of the whole structure. This requires a full 
derivation of the multilayer structure transmission in the 
TF configuration, a rather complicated problem which 
we have succeeded to solve in some specific cases and 
which will be detailed in a further publication ^(|. We 
just here outline the principles. We used a standard ap- 
proach which consists in evaluating for each layer N, with 
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a dielectric tensor e n(uj), the transfer matrix Mm after 

finding the appropriate modes of propagation of light k N 
inside the layer, by solving the Maxwell equations. Tak- 
ing the z-axis perpendicular to the layer, and defining the 
y-z plane as the plane of incidence which also contains 
1? , the angle 9 is defined as the angle There 
are four modes k\ N for each layer. The 4x4 transfer 
matrix is built by writing the conservation of the tangen- 
tial components of the electric and magnetic fields of the 
light at each interface of the layer. If e ^(uS) is diagonal, 
like in layers which do not contain free electrons (or at 
B = OT) the transfer matrix decomposes into two blocks 
of 2 x 2 matrices related to the TE and TM modes re- 
spectively. This result is well known and has been used 
to measure the LO phonons frequencies of semicon- 
ductors in a way similar to the results reproduced in Fig. 
la. For the doped QW, in the TF configuration, all el- 
ements of e qw(u) and therefore of the transfer matrix 
are non zero and the TE and TM modes are mixed. The 
resulting total transfer matrix has the same properties 
and a special treatment 0] needs to be developed to 
extract the transmission of the structure. This can be 
done as soon as the dielectric tensor of the QW in the 
TF configuration has been defined. 

In the PF configuration, for any finite value of B, the 
z-part of the 2DEG wave function is decoupled from the 
x-y part. When 9 / 0, this decoupling is no longer valid 
and as we shall see below is indeed responsible for the 
observed anticrossing. One has then to reproduce the 
mixing of the z and x-y part of the 2DEG wavefunction 
with the best accuracy. For the square QW which is un- 
der study, this mixing can only be evaluated in perturba- 
tion theory On the other hand, if we assume the z- 
confining potential to be parabolic, with a characteristic 
energy TiQ, separating the electric subbands, the coupling 
can be evaluated analytically at all orders. Of course the 
strength of the coupling is modified, the dipole oscilla- 
tor strengths along the z-axis being different between a 
square QW and a parabolic one. It is this model, based 
on the Drude formalism, that we have used to evaluate, 
neglecting retardation effects, all components of £qw(u) 
and then simulate the transmission. Before discussing 
these simulations it is instructive to look at the struc- 
ture of two characteristic components of £ qw{oj)- The 
real part of e rx and e zz are expressed, ignoring damping 
effects, as TO] : 
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FIG. 3: Fit of the relative transmission curves for sample 
1200 at 8 — 25°. The curves are shifted for clarity. Empty 
dots are experimental points and full continuous curves the 
corresponding calculated relative transmission spectra. The 
hatched region corresponds to that of Tiuito (GaAs) (see text). 

and u>2 ar e the solutions of the equation: 



w 4 - (ir + w^y + tui ± n z = o 



(3) 



where El is the contribution of the lattice, <Ji = 
4ire 2 ns/(Lm*) the square of the plasma frequency, ui\ 



with lu c = eB / (Tim*). When Q >> uj c ±, oj\ — > u c ±_ 
and 102 - * ^- In the present case, however, the exact 
values for the poles are introduced in the calculation. 
This explains why in Fig. 2, far from the anticrossing 
region the energies of the CR lines do not coincide with 
the values obtained in the PF configuration. 

Different fitting steps have to be performed. We first 
have to fit El with standard values as already done 
and shown in Fig. la (red continuous lines). In fact, in 
the region of Tiluto (GaAs), interference effects, between 
phonon and CR absorptions become important and this 
is the reason why we prefer to interpret the relative trans- 
mission spectra (Fig. lc and Fig. 3) for which these spu- 
rious effects are minimized. They are however not elimi- 
nated and, for this reason, we ignore in the fitting process 
a region of ±0.5 meV around Hcuto (GaAs), sketched as 
the hatched region in Fig. 3. In this figure, experimental 
data (empty dots) and their simulations (red continuous 
curves) for sample 1200 are displayed for 9 = 25°. There 
are different fitting parameters to adjust: ns, Tri*, tqb, 
, O and 9, tqr being the scattering time entering the 
Drude model . In the PF configuration, rig is fitted at 
lower fields with the model extended for NP effects [lcj 
together with m* and tcr- Outside the energy region 
around Tiloto (GaAs) where there are signs of some in- 
teraction (at least for the sample 1200) which are not 
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discussed here, the values of m* fitted correspond quite 
well to the expected variation induced by NP effects [l3| 
whereas the mobility )j,cr = eT CR.I m * ranges around 
150 m 2 /V/ sec revealing the good quality of the samples. 
The angle 6 has also to be fitted because, for technical 
reasons, the Si substrate being wedged, 8 can differ by 
a few degrees from the mechanical angle. It is fitted by 
setting the m* value constant for a given value of B±. 
This is not completely independent of the value of fl but 
the loop easily converges. The value of B± indicated in 
the figures are calculated values of B x cos{0). Finally 
il is fitted in such a way the splitting of the CR1-CR2 
lines can be best reproduced. We impose also, in the fit- 
ting process, the value of to be constant for all angles 
and QW having the same width L. As an example, the 
results, shown in Fig. 3 for sample 1200, are obtained 
for a value of f2 = 63 meV (with an uncertainty of ±3 
meV). These results clearly demonstrate the anticrossing 
features observed experimentally. The fit reveals some 
slight deficiency which is sample dependent but the over- 
all agreement is surprisingly good. It is interesting to 
note that if one calculates, at B = 0T, the zeros of e zz 
(Eq. 2) with the fitted value of fl, one finds for the low 
energy modes 35.62 meV for the sample 1200 and 35.49 
meV for the sample 1211, values which agree very well 
with the experimental findings. These modes are the so 
called plasmon-phonon-intersubband (PPI) modes (lif . 
Therefore the observed anticrossing occurs between the 
CR mode and this hybrid PPI mode. They are coupled 
by symmetry but not by any specific electron-phonon 
interaction. We think, though approximate, the model 
keeps all the physical and symmetry aspects of the prob- 
lem. One direct consequence is that, when increases, 
the oscillator strength of the PPI-like mode rapidly de- 
creases to zero ll(J. Therefore, for a pure 2DEG, there is 
no longer any interaction whatever is the angle. The data 
in the PF and TF configurations do not show any sign 
of interaction related to the Frohlich coupling: we are 
therefore led to conclude that the concept of polaronic 
mass, related to this interaction, is no longer effective. 

All the discussion is based on the assumption that we 
are dealing with carriers for which a plasma frequency 
can be defined. It is known that for bound electrons, 
like the ones found in neutral shallow impurities, the 
Is — 2p + absorption reveals 0] the Frohlich interaction 
around the energy of the LO phonon in GaAs but then 
the mass of such an electron is not defined. The same 
situation holds for electrons in quantum dots It is 
also true that there is still some electron-phonon inter- 
action present in such system : indeed, around the TO 



frequency, there is a clear sign of interaction especially in 
sample 1200. This is an interaction which may be driven 
by the mechanism of deformation potential which is at 
the origin, for instance, of the Raman response of such 
structures. 

In conclusion, we have performed careful magneto- 
optical measurements of the infra-red absorption related 
to quasi two-dimensional electron gas, with different den- 
sities, in a range of magnetic fields which allows the cy- 
clotron resonance to span the optical range of phonon en- 
ergies. All singularities observed can be explained quan- 
titatively without any electron-phonon coupling of the 
Frohlich type. We argue that this support the idea that 
the concept of the Frohlich polaron mass has to be re- 
examined in a real material. 
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